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D ■ Abstract 

. To date, the properties of Q-balls arising from an Affleck-Dine condensate in 
gravity-mediated SUSY breaking have been obtained primarily through numerical 

^ ' simulations. In this work, we will derive the expected charge of the Q-balls formed 

. in such a scenario through an analytical treatment. We will also examine the numer- 

\ ically observed difference between Q-ball formation in weakly charged condensates 

^SJ ' and formation in strongly charged condensates. 
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1 Introduction 



D , 

■ A Q-ball is a non-perturbative solution of the equation of motion for a scalar 

field which is charged under a continuous U(l) symmetry [1]. The MSSM 
^ . requires several such scalars, and so it is expected that Q-balls could be formed 

in a supersymmetric universe [2]. An obvious context in which to look for Q- 
balls is the (baryon charged) scalar condensate that is required for Affleck-Dine 
(AD) baryogenesis [3]. This possibility has been explored both analytically 
using linear perturbation theory and numerically. 

On the analytical side, two distinct methods of formation have been discussed 
in the literature. A strongly charged condensate might fracture into Q-balls. 
This possibility is discussed in [4]. It has also been noted that the character- 
istic AD scalar potential in gravity mediated SUSY breaking can give rise to 
negative pressure that causes perturbation growth [5]. 

Numerically, exhaustive simulations have been performed by Kawasaki and 
Kasuya [7,8]. They observe Q-ball formation for a wide range of initial AD 
field configurations, including both strongly and weakly charged condensates 
(we will quantify the terms "strongly" and "weakly" later in the paper). These 
simulations reveal that the mechanism of Q-ball formation is not the same in 
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all condensates. For strongly charged condensates, a proportionality between 
the maximum Q-ball charge and the condensate charge density is observed. 
There is, however, a distinct change in this relationship for weakly charged 
condensates. This change in scaling behavior is accompanied by (or results 
from) the appearance of Q-balls carrying a charge opposite in sign to that 
initially carried by the condensate in addition to those of the same sign. Thus, 
for weakly charged condensates, numerical simulations imply that Q-balls can 
be effectively pair-produced. 

In this work, we will briefly revisit the analytical work that has been done 
on Q-ball formation. We will show that the two different approaches to Q- 
ball formation can help to explain why weakly charged and strongly charged 
condensates form Q-balls differently. We will attempt to explain the observed 
relationship between the charge density of the condensate and Q-ball charge 
for strongly charged condensates. 



2 Review of Q-balls in the AD Scenctrio 

2.1 Flat Directions 

AD baryogenesis will require a scalar field that carries baryon number and has 
a very flat potential so that it can be given a large expectation value [3] . There 
arc numerous flat directions in the renormalizable MSSM without soft SUSY 
breaking terms. These have been cataloged [9]. High energy (> Mgut) physics 
and SUSY breaking, however, ensure that no direction will be completely flat. 
In modern formulations of the AD scenario [10,11] non-renormalizable, soft- 
breaking, and inflaton-coupling terms are added to the renormahzable MSSM 
potential. Such terms lift a flat direction by generating a potential of the 
form: 

U{<t>) = K - ^hH'M' + ^.M"- + h.c.) + ^^101^"-^ (1) 

where cj/ ~ 1, |am"^3/2| ~ fn^ ~ 1 TeV (assuming gravity-mediated breaking), 
and M is a large mass scale. 

The phase dependent term is necessary for AD baryogenesis, as it will allow 
dynamical charge creation in the condensate. For the purpose of this dis- 
cussion, however, it is irrelevant. We will simply assume charge creation has 
occurred. We are instead most interested in the negative mass term propor- 
tional to the Hubble constant. If we have a large enough H during inflation, 
the (f) fleld will spontaneously move out to a flnite expectation value. (Note 
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that the negative sign is chosen. It is equally likely to be positive [10], but this 
is an uninteresting case.) 



— 1/2 

Once H drops below the value Hgsc — 'm^Cfj , however, we expect the positive 
mass term to begin to dominate, and will begin coherent oscillations about 
zero. We can find the value of when the oscillations begin. For H ~ i^osc, 
minimizing the potential (1) yields a (j) vacuum value of: 

, ^ l/(n-2) 

|A|2(n-l)y 



This gives a characteristic initial condition for the scalar condensate. Q-ball 
formation, however, has to wait until the phase- dependent term has become 
irrelevant, otherwise charge is not conserved. This will require the field value 
to drop by a factor of order 10 from the value of equation (2). Examining the 
form of the potential (1) tells us that for such small (f) values, we can use the 
effective baryon-conserving potential: 

U{cj>)=ml{m? (3) 



where the 4> dependence of is essential to Q-ball stability. We will discuss 
this point in the next section. 



2.2 Radiative Corrections and the Q-ball Condition 



The next step in forming Q-balls is to show that the potential (1) is compatible 
with the Q-ball condition [1]: 



mm 
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<-u"m = o). 



(4) 



This condition can be realized when radiative corrections are taken into ac- 
count [5] . The gauge interactions associated with the components of the (com- 
posite) field add a scale-dependence to the mass, so that: 



^ mlifj,) l + Xlog 



(5) 



where K is negative, with a magnitude in the range 0.01 to 0.1 [5,6]. The 
negative value of K is important, as it implies that equation (4) is satisfied. 
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The potential relevant for Q-ball formation, then, is: 



[/(|0|) = mJ(/.)|0|Ml + i^log 



+ 



lAI 



2n-2 



(6) 



Given the form of this potential, the solution for thin-walled Q-balls has a 
characteristic (f) value [1,5]: 
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(7) 



Note that this is almost the same expression as our expected initial value of 
(f) when the condensate begins oscillations (equation (2)). This is something 
of a problem for Q-ball formation. Recall that Q-balls cannot form until the 
phase dependent term has run its course. We do not expect the AD field to 
be large enough at that time to form thin-walled Q-balls. 



2.3 Thick- Walled Q-balls 



Since AD baryogenesis will have to form thick-walled Q-balls if any, we must 
consider the properties of such constructions. These were examined in [6]. 
There it was shown that we can assume a Gaussian profile for thick walled 
Q-balls. That is, we can take: 

(j){r, t) = 0(r)e^'"* = |0e|e-'''/^'e*'^* (8) 

where 4>c and R are constants describing the radial profile of our Q-ball. 

To see why this works, we must consider the equation of motion for the radial 
part of (f) [1] : 



gj.2 J. Qj. Q^* 



(9) 



For a radial profile of the form (8) and a potential of the form (6) this yields: 



6 2 2 ir^l 2 



+ (10) 



where we have ignored the contribution of the non-renormalizable term on the 
right hand side, because for 0c < 0o it will be negligible. 
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Comparing the sides of the above expression term by term then tells us: 



\K\ml 



while 



uT ^mi + 2\K\mi 



;i2) 



We can also determine the charge of a Gaussian Q-ball. The Noether current 
for our complex scalar is: 



(13) 



where f3 is the baryon charge per meson of the field 0. Now we can use equa- 
tions (8) and (13) to see that the total charge can be expressed [6]: 



3/2 



(14) 



3 Strongly Charged Condensates 



3.1 Background 



Having studied the properties of Q-balls, we now investigate whether they are 
likely to be formed in an AD baryogenesis scenario. For a strongly charged 
condensate, we anticipate that the uniform scalar field could simply fracture 
into stable Q-balls. In the last section, we showed that the Q-balls will have 
a characteristic size. Therefore, the process of breakup should yield Q-balls 
with a charge Q proportional to the initial charge density of the condensate 
go [4] where qq in the simulations of [7,8] is defined as the charge density when 
H pa rrifj). 

The exact proportionality will depend on the time at which breakup occurs. 
This is because in an expanding universe the baryon charge density will fall 
with time. For a matter-dominated universe (Affleck-Dine baryogenesis occurs 
during the epoch of inflationary reheating when inflaton matter dominates the 
energy density of the universe), the charge density decreases as g oc oc H^. 
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3.2 Q-ball Stability 



We first show that a strongly charged condensate is energetically capable of 
breaking into stable Q-balls. Recall that Coleman's condition for Q-ball sta- 
bility is that the energy per charge in the Q-ball be less than the energy 
per charge in free particles of the scalar field in question. The energy of free 
scalars, however, is simply the mass of the field: = U"{0)/2. This line of 
reasoning leads directly to equation (4) [1]. That equation tells us that the 
binding energy per charge for a thin-walled Q-ball will be: 




where 0c is the field value within the Q-ball. Note that this binding energy can 
be maximized in gravity mediated scenarios at the specific value 0o defined in 
equation (7). 

A problem arises when we try to calculate the value of the binding energy, 
however. It is a well-known difficulty (see e.g. [12]) that the running mass of 
our scalar field has an infrared divergence associated with = 0. Even worse, 
however, the one-loop running becomes untrustworthy for large values of the 
log term, since the next loop order will bring in higher powers of the log. This 
implies that we cannot actually calculate the mass of a free meson by one- 
loop running from the value of the mass at large and, more importantly, the 
qualitative behavior becomes suspect at large values of the log term. 

This could be a very important issue for thick-walled Q-balls, which necessitate 
0c << 00 SO that we are already stretching the one-loop running to the limits 
of applicability. If we want to have a truer idea of the Q-ball binding energy, 
we should try to make use of the renormalization group, or at the very least 
move to two-loop running. This has not been attempted in the numerical 
simulations to date, however. Therefore, let us consider the consequences of 
assuming simple one-loop running holds over a large range of 0. 

With this assumption, we see that although 0o maximizes the binding energy 
per charge, the binding energy will remain positive (formally infinite) for any 
value < 00. This is very important for us, since we have already shown 
that the AD scenario will only generate thick-walled Q-balls. We can now see 
that they are potentially stable. Unfortunately, our discussion so far is not 
technically valid for thick-walled Q-balls since they will require the addition 
of a gradient term to the expression for the Q-ball energy. 



6 



The energy of a complex scalar field is given by [1] : 



E 



(fx 



+ V0* + t/(|0|) 



(16) 



Note that if we do have gradient energy, we can treat it as an effective con- 
tribution to the scalar potential. That is, we can define an effective potential 
for thick-walled Q-balls: 



= + V0 • V0* 



(17) 



We can estimate the contribution of the gradient term to U by using the 
Gaussian approximation. Recall that thick-walled Q-balls have the profile: 



</.(r,t) = |0,|e-i"l"^^ 



(18) 



which means: 

V(/) • = |X|2m>2e-I^K'-'|0/. (19) 



Looking at this equation, we see that the gradient term will make a contri- 
bution to the potential which is at maximum of order |iir|m^|0cP- Since the 
energy per charge has the form: 



E 



Q P\ 



(20) 



we expect the gradient energy to yield a contribution of order \K\^/'^m^. We 
would like to know if this is significant on the scale of the expected binding 
energy per charge. 

Even though the binding energy is formally infinite, we can still compare the 
characteristic scale of the binding energy with the characteristic scale of the 
gradient energy. In the thick- walled case, we expect that the potential will be 
dominated by the mass term. Thus, we can estimate the binding energy per 
baryon as: 

^y^- ^V^^(0c) ^ ^|K|i/2m<^log(0e/O). (21) 

This equation shows that although the binding energy is formally infinite, 
it has the characteristic magnitude \K\^/'^m^ (we will assume that the log 
running is somehow regulated by higher-order effects). This implies that we 
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can expect a binding energy per charge of at least a few times \K\^^'^m^, 
and so the gradient energy of thick-walled Q-balls is not too large to admit 
bound Q-balls. These conclusions are borne out in numerical simulations. In 
[7,8] Q-balls are shown to form over a wide range of initial values (many 
decades) . 



3.3 Partially Charged Q-balls 

So far we have only considered Q-balls which follow a circular path in phase 
space (a; is constant). As we shall see, this corresponds to a scalar field that is 
completely charge asymmetric (all baryons or all anti-baryons). This is not a 
necessary condition for a strongly charged condensate, however. In [7] it is seen 
that even condensates with considerable ellipticity in their phase space rota- 
tion will form Q-balls that obey the same scaling relations as for completely 
charged condensates. This implies that elliptically rotating condensates can 
fracture into bound Q-balls as well. We would like to examine the energetics 
of this process. 

The first step is to parameterize an elliptical condensate. Recall our effective 
potential after the phase dependent terms have become unimportant is: 

C/(|0)=m,|0r(l-|ir|log(^|)) (22) 

where 0i is basically arbitrary. As long as does not change more than an 
order of magnitude or so, the log corrections are essentially unimportant as- 
suming small \K\ (as we shall see, this is true for a condensate with significant 
charge). Therefore we effectively have a |0p potential. Such a potential allows 
for closed orbits in the complex plane [15]. We can parameterize such an 
orbit in the form: 

(f) — Acos{m^t) + iSsin(m^t). (23) 

Note that in keeping with our orbit analogy the analog of angular momen- 
tum in this case is which, by the Noether current, is equal to nB/2/3 
where is the number density of excess baryons. Conservation of angular 
momentum is guaranteed by the central force character of our potential (even 
with logarithmic correction). Plugging our parameterization into the Noether 
current yields: 

riB = 2l3m^AB. (24) 
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We can now compare this number to the total number of particles that our 
condensate contains. This is approximated by taking the total energy density 
in the condensate, as given by equation (16), and dividing by m^: 

m^A^ + m^B'^. (25) 



It is now apparent that if we have a circular "orbit" the condensate is entirely 
composed of baryons (or antibaryons) . For elliptical orbits, higher eccentricity 
means less net baryon fraction. 

Now let us evaluate the energy per charge of an elliptical Q-ball. We can 
make the analysis easier by picking a particular point along the trajectory. 
The simplest to use is when (f) = A. Using our small (f) approximation for the 
potential and ignoring gradient terms for simplicity, we obtain: 

E ^ {mlB^ + mlA^)V (26) 



where V is the volume occupied by the Q-ball. The charge of the Q-ball will be 
ubV where is given in equation (24). Thus, we have an energy per charge: 

E I f B A\ 



Note that throughout this analysis we have ignored the running of the mass 
with (f). Now we can see that this was justified. The characteristic contribution 
of the ellipticity to the energy per charge is of order m^, while the binding 
energy is characteristically of order \K\^/'^m^. 

This is a serious concern. By comparing equation (27) to the result of equation 
(21) we can see that once A is a few times greater than B it would require 
the one-loop mass running to be vahd over several orders of magnitude to 
yield a bound Q-ball. In contrast, numerical simulations show that elliptical 
condensates can certainly form Q-balls up to ratios A/ B ~ 10, and possibly 
beyond [7]. 

We can reconcile this apparent discrepancy by noting that the cause of the 
instability of Q-balls with a large energy per net baryon is quantum mechan- 
ical decay. Such decay is not treated in numerical simulations of the classical 
field. Thus, our estimates would lead us to conclude that unless some way of 
damping out the ellipticity of a partially-charged AD condensate is present, 
we must question whether such condensates will truly form Q-balls. 

One possibility for such damping is scalar self-interaction in a nonhomogeneous 
condensate. It is possible that this has already been observed in the numerical 
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work of [7,8] and simply not published. Perhaps a more interesting possibility, 
however, is that in this case the problem is also the cure. 

A rotating scalar condensate need not decay only to free scalar particles. We 
anticipate couplings of the AD field to many other particle species. In this 

case, we would anticipate that a dense condensate could experience annihi- 
lations that remove the charge-neutral part while leaving behind the baryon 
asymmetry. This would effectively remove the excess energy per baryon. 

We have calculated the annihilations that would be mediated by an interaction 
term of the form where x is another scalar field. Such an interaction 

would generically be present for a flat direction composed of squarks. Our 
calculations show that the annihilations mediated by this interaction would 
reduce the ellipticity of the condensate dramatically on a time scale that is 
considerably less than the time scale for Q-ball formation [16]. 

3.4 Breakup 

So far in this section, we have summarized the evidence that a uniformly 
charged AD condensate is in a state that meets the criteria to be considered a 
bound Q-ball immediately after beginning coherent oscillations about = 0. 
On the strength of these arguments, we expect that by the time the universe 
is of the age t ~ 0{1) x rri^^ the condensate is energetically ready to form Q- 
balls. Numerical simulations, however, show that the time of Q-ball formation 
is more nearly of the order t ~ 0(1000) x m^^ for \K\ = 0.01 [7]. 

This discrepancy is a result of the finite age of the universe (or, equivalently, 
the expansion of the universe). In a matter dominated universe, the maximum 
radius of a causally connected patch is 

r,, = 3i = 2H-\ (28) 

Thus, even if the scalar field would energetically prefer to form Q-balls, the 
Q-ball sized patches cannot form a coherent object until Vcc ~ \K\~^^'^m'^^. 

For \K\ = 0.01 Q-ball formation would then be causally allowed by t ~ 
C(10) X m^^. Thus, we still have a large factor to make up. As we shall 
see, this is because the forces responsible for Q-ball formation propagate at a 
characteristic sound speed rather than the speed of light. 

Q-ball formation is the result of interactions in the fluid of squarks represented 
by the scalar condensate. Such fluid interactions will have a sound speed asso- 
ciated with them. The sound speed in a completely charge-asymmetric scalar 
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condensate is given in [1] as: 



2 ^ u"m\<i>\-u'{^.u ..Qi 



Evaluating this using the effective potential of equation (3) gives the result 
that v"^ K/2. It is important to note that because K is negative, Vg is 
technically imaginary. As we shall see, this will lead to the amplification of 
small perturbations. 

To see this amplification, consider the plane- wave solution to the wave equa- 
tion: 

ijj{f,t) = ijjoexp (ik ■ f±i(jJst) (30) 



where for sound waves {u;s/\k\y — v^. 

Now, consider the situation when is negative. We can construct a very 
similar solution, except that now we must take cus to be imaginary. Thus, we 
anticipate solutions of the form: 

-0(7^, t) = ipQ exp iik ■ r \ exp (ila;^!^) . (31) 



This solution yields growing modes which increase exponentially with e- folding 
time equal to: 

I 1-1 1 

te-foU-\u^s\ (32) 



Since a sound wave is constructed of small density fluctuations in the medium, 
this tells us that such fluctuations will grow exponentially. Further, since the 
density of the medium is proportional to we anticipate that small fluc- 
tuations in the field behave as 5p ~ so that fiuctuations will grow on 
basically the same time scale as density perturbations. 

Consider what this means for our condensate. We expect to form Q-balls with 
characteristic diameter 2^/^|i^r|~^/^m^^. Thus, to decide how fast perturba- 
tions on this scale will grow, we can use this diameter as half the wavelength 
in determining \k\va. the above equation. We have already calculated Vs- Com- 
bining these results with the results above gives: 

2%/2 1 , , 

te-fold ~ -T7T\ ~ 7-rr\ • 33) 
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For \K\ = 0.01, this evaluates to a characteristic time scale for the growth of 
Q-ball-sized perturbations: 

te-foid ~ lOOm.^ (34) 



Further, since the numerical simulations of interest start off with perturba- 
tions suppressed by a factor of 10^ [7], we anticipate Q-ball formation will 
require about 16 e-folding times. This yields exactly the relation found in the 
numerical simulations: 

tq-haii ~ iQte-foid ~ SOOOm^^ (35) 



Calculating the total charge contained in the Q-ball sized perturbation at this 
time, we find that: 

Q ^ -jRy{0)tQl,,um/ « 0.0057m^3^(0) (36) 



which compares very well with the results obtained on three-dimensional lat- 
tices in [8]. 

Before we leave this section, we should comment that this analysis is only 
valid for imaginary values of the sound speed in Q-matter. Looking back at 
equation (29), this means that Q-ball formation will require: 

u"m-u'mm<o. (37) 



This condition can be restated using the equation of motion for From this, 
we can see that the phase rotation rate cu for a completely charge-asymmetric 
condensate is given by [1]: 

u^'^U'm)M. (38) 



Thus, we can rewrite our condition for Q-balls: 

U"{\(P\)-u'' <0 (39) 



which is precisely the condition for Q-ball formation in a completely charge- 
asymmetric condensate arrived at in [4] using the equations of motion for |0| 
and cu. 
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4 Weeikly Chcirged Condensates 

4.1 Background 

For condensates with little or no net charge, we expect that breakup into Q- 
balls must follow a different path than for strongly charged condensates. The 
only way to form Q-balls is to somehow separate the neutral condensate into 
positive and negatively charged regions. This is observed to happen in numer- 
ical simulations [7]. We wish to try to understand this process analytically. 

4-2 Negative Pressure and Growth of Perturbations 

In one of the early papers on Q-balls in gravity-mediated AD scenarios, it 
was assumed that Q-ball formation was due to negative pressure in the AD 
condensate [5]. For the case of a real (non-complex) scalar field, negative pres- 
sure condensates result from the negative value of K in the one-loop corrected 
scalar potential (see equation (6)) [13,14]. We wish to exphcitly consider the 
pressure arising for a charged (complex) scalar. 

It is important to note that in this section we will have to be careful to retain 
terms of order \K\ in the expressions we are about to derive. These small 
corrections will be responsible for any pressure that we find. The first such 
correction to take note of is that when we parameterize our rotating condensate 
in the form of an ellipse, we must use the more general form: 

(f)^ A cos{ut) + iB sn\{ijjt) (40) 
where to first order in a; is given by: 



(This result follows from the equation of motion.) Note that for a treatment 
of "weakly charged" condensates we are interested m. A/ B > 10. 

With this form in hand, we can move on to explicitly calculate the pressure. 
In [13] it was shown that for a rapidly (relative to the Hubble time) oscillating 
condensate, the average pressure can be written in terms of the energy density 
through the use of the parameter: 




(41) 




(42) 
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where p is the pressure, p is the energy density, and T is the period of the 
oscillation. In terms of 7, we would have the relationship between average 
pressure and energy density: 

p=(7-l)p. (43) 



To evaluate this expression for our charged condensate we must use the energy- 
momentum tensor for a complex scalar: 

^'''=m''^^m-'''^'-''''' '''' 



which tells us: 

p = + [/ (45) 

and: 

p = 00* - U. (46) 

These, along with the generic elliptical parameterization of equation (40), yield 
the expressions: 

p + p = 2#* = 2(p-C/) (47) 

and: 

p = ml [A" + B^) - \K\mlB\ (48) 

Now, assuming that p changes very little over the course of one oscillation, we 
can write: 



2-k/uj 

^ 



This expression can be numerically integrated. Doing so for various values of 
B/A and K yields the curve shown in Figure 1. The feature of current interest 
is the intercept at S = 0, which is consistently equal to — 0.39|ir|. 

With this result we see that the analysis of [5] is essentially correct. The 
pressure of a weakly charged AD condensate is expected to be negative, with 
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(7-1)/|K| 



-0.55 



-0.45 



-0.4 
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Fig. 1. Plot of p/{\K\p) for condensates with varying ellipticities. The sohd hne is 
K = -0.1, while the dotted line is K = -0.01. 

7 — 1 ^ K/2. Note, however, that we have made something of an approx- 
imation in our parameterization of the AD condensate as a perfect ellipse. 
A study of "true" AD condensates obtained via numerical integration of the 
classical equations of motion shows that the numerical treatment for the ide- 
alized elliptical trajectory is accurate to within a factor of 2. Therefore, we 
have settled on K/2 as a reasonable first approximation to the pressure. Note 
also that this is the result obtained using the polynomial approximation of 



We can now follow [5] to see that perturbations in even a perfectly charge- 
symmetric AD condensate will grow exponentially. In fact, since the factor 
7— 1 is identical to the sound speed Vg found for a strongly charged condensate, 
the exponential growth will occur with the same timescale as the growth of 
perturbations in a strongly charge-asymmetric condensate. 

4.-3 Charge Separation 

The negative pressure of the condensate can serve to give us growing Q-ball 
sized perturbations. The next question, however, is how these regions can 
acquire a charge asymmetry. The dynamics of such a process are beyond a 
simple perturbative analysis. We can, however, do some back of the envelope 
calculations to guess at the charges we could expect. 

Consider two overdense Q-ball-sized perturbations that have evolved "next 
to" one another. We will consider the epoch when these perturbations have 



[14]. 
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just gone nonlinear, so S = Sp/p ~ 1 where p is the average energy density 
in the scalar field and Sp is the peak energy density in the perturbation. The 
perturbations will be separated by an underdense region of approximately 
Q-ball size. 

We must assume that they have interacted in such a way as to impart on one 
another a very weak rotation. By conservation of baryon number we know one 
will have a rotation corresponding to positive baryon asymmetry and the other 
will carry negative baryon asymmetry. We now wish to show that this weak 
rotation is enough to cause the positive perturbation to preferentially accrete 
baryons as it grows, while the negative perturbation will preferentially ac- 
crete anti-baryons. Thus the charge asymmetries will be enhanced, eventually 
yielding a Q-ball/anti- Q-ball pair. 

Consider a tiny box (side length L) filled with baryons that lies on the line 
joining the two perturbations. Imagine that it is closer to the center of the 
negative perturbation. We can calculate work required to remove this box from 
the negative perturbation. The work will be: 

oo 

Wrm = j ApL^dr (50) 



where ro is the distance between the center of the negative perturbation and 
the small box to begin with, and Ap is the pressure difference between the 
leading edge and the trailing edge of the box, given approximately by 

Ap ^ -^L. (51) 
dr 

So, using the results of the previous section: 

ro 



We want to compare this work with the energy difference that we would expect 
to see if the box of baryons left the negative perturbation and joined the 
positive perturbation. We can estimate this by considering the perturbations 
to be highly elliptical Q-balls. Thus, we expect that the perturbations contain 
net baryon number Nb equal to (see equation (25)): 

Nb = ±2Pm^ABVQ (53) 



where the minus denotes the negative perturbation and Vq is the volume of 
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the perturbation. Similarly, each perturbation has an energy (see equation 
(26)): 

E = {mlA^ + mlB^)VQ (54) 



where in each of these formulas, we expect B. 

Adding the box of baryons to either perturbation implies a change in both the 
total baryon asymmetry and the energy. For the net baryon number, we have: 

5Nb = /3— = ±2/3m^ {A5B + B5A) Vq (55) 



where pi, is the energy density within the small box, and we have differentiated 
Nb given above. The energy relation is: 

5E = pi,L^ = 2m! {A5A + B5B) Vq. (56) 



We can solve these equations simultaneously to obtain: 
AR Pfci^' A-B 



for the positive perturbation, while for the negative perturbation we would 
have: 



We have solved for these by explicitly assuming the energy contribution of 
the box of baryons is the same for each perturbation. To zeroth order in K 
we expect this to be the case, as there is no attractive force present. Once we 
"turn on" K, however, we do expect to see an attraction between baryons. We 
can estimate this attraction by using our zeroth order solutions for SA and 
SB in the first order expression for energy. From equation (48) we can see: 

BA B 

EunUke - Euke = 2\K\p,L^ ^ 2\K\p,L' -. (59) 



where "like" denotes the energy of a positive perturbation with the positive 
box of baryons or of a negative perturbation with a negative box, and "unlike" 
represents the other two possibilities. Thus we see that at first order in K the 
energy of the box is lower when it is within the positive Q-ball than when it is 
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within the negative Q-baU. This energy difference can provide the work needed 
to separate the negative charges from the growing positive perturbation and 
vice- versa. 



4-4 Q-ball Charge 

We will now make a very rough estimate of the amount of charge separation 
that will occur in the formation of a perturbation from a weakly charged 
condensate. 

In the simplest picture, we are forming spherical perturbations with radius R 
(approximately |i^|~^/^m^^) which are separated from one another by voids 
of about 2R in width. Thus, each perturbation has a sphere of influence that 
extends a radius 2R from its center. Theoretically it will gather up all the 
matter in this sphere to create the final perturbation. 

We will now assume that we have a pair of oppositely charged perturbations 
next to one another. Now in the region between these, the baryons will be 
preferentially attracted to the positive perturbation, while the antibaryons 
will be preferentially attracted to the negative one. Thus, we will somewhat 
arbitrarily define an extended sphere of infiuence for these charged Q-balls 
of radius SR. In the overlap region of these extended spheres, we will assume 
that the charges separate completely and collect in their energetically preferred 
Q-ball. 

To determine the final charge, then, we are interested in the ratio of this over- 
lap volume to that of the complete 2R sphere of infiuence. A quick calculation 
shows that the ratio in question is approximately 1/10. Thus, we expect a 
final charge fraction in the Q-ball of 1/10 the value that a Q-ball formed in a 
completely charge asymmetric condensate of the same initial |0| value would 
have. This is in reasonable agreement with numerical simulation [7,8]. 



5 Summciry 

We anticipate that strongly charged condensates can easily break up into 
bound Q-balls. It is clear from a simple analysis, however, that any significant 
ellipticity in the phase rotation of the condensate might compromise this bind- 
ing by making the energy per charge larger than that in a free scalar plasma. 
This concern leads us to believe that damping of this ellipticity is important 
for Q-ball formation. This has been considered in a separate work [16]. (Note 
added: After publication of this paper, I was informed that the delayed decay 
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of Affleck-Dine condensates has also been considered in [17].) 

Assuming that strongly charged condensates can reach a nearly circular ro- 
tational state, we recover the known result [4] that breakup will occur. By 
considering the sound speed in the scalar condensate, we have shown that 
the time scales for breakup observed in numerical simulations [7,8] are repro- 
duced in a simple analytical treatment. Given the results of [7,8] we expect 
these "strongly charged" results to be useful for condensates with A/B a.s high 
as 0(10). 

For condensates with even larger A/B ratios, we have shown that the neg- 
ative pressure of the condensate [5,14] can still allow for the growth of seed 
perturbations. This, coupled with the weak attractive force between baryons 
can lead to Q-ball/anti-Q-ball pairs. 
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